Control of spin relaxation in semiconductor double quantum dots 
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We propose a scheme to manipulate the spin relaxation in vertically coupled semiconductor double 
quantum dots. Up to twelve orders of magnitude variation of the spin relaxation time can be achieved 
by a small gate voltage applied vertically on the double dot. Different effects such as the dot size, 
barrier height, inter-dot distance, and magnetic field on the spin relaxation are investigated in detail. 
The condition to achieve a large variation is discussed. 

PACS numbers: 73.21.La,71.70.Ej,72.25.Rb 
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Spin related phenomena in semiconductor nanostruc- 
tures have attracted much interest recently due to the 
fast growing field of spintronics 1 . Among different struc- 
tures, quantum dots (QDs) have caused a lot of attention 
as they provide a versatile system to manipulate the spin 
and/or electronic states 2 . Many proposals of spin qubits, 
spin filters, spin pumps and spin quantum gates are pro- 
posed and/or demonstrated based on different kinds of 
QD S 2, 3,4,5, 6,7,8,9,10, ii _ Manipulation and understanding 

of the spin coherence in QDs are of great importance in 
the design and the operation of these spin devices. There 
are many theoretical and experimental investigations on 
the spin relaxation in single QDs 12 ' 13 - 1445 - 16 ' 17 - 18 , double 
QDs 19,20 and quasi-one-dimension coupled QDs 33 due to 
the Dresselhaus or Rashba spin-orbit couplings 22,23 . In 
this paper, we propose a feasible and convenient way to 
manipulate the spin relaxation in double QDs by a small 
gate voltage. We show that up to twelve orders of mag- 
nitude variation of the longitudinal spin relaxation time 
(SRT) can be tuned in such a system. 

We consider a single electron spin in two vertically cou- 
pled QDs. Each QD is confined by a parabolic potential 
V c (r) = |m*Wjr 2 (Therefore the effective dot diameter 
do = \/ fnr / m* loq) along the x-y plane in a quantum well 
of width d with its growth direction along the z-axis. 
A gate voltage Va together with a magnetic field B are 
applied along the growth direction. A schematic of the 
potential of the coupled quantum wells is plotted in the 
inset of Fig. 1(a) and the potential is given by 24 
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in which Vo represents the barrier height between the two 
coupled QDs, a is the barrier width and E = Vd/(a + 2d) 
denotes the electric field due to the gate voltage. The 
origin of the z-axis is chosen to be the center of the bar- 
rier between the two QDs. By solving the Schrodinger 

equations along the z-axis d 2 ip z / d£% — £,iip z = with £i = 
2 i/3 ( 
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FIG. 1: (Color online) (a) SRT vs. the electric field. Solid 
curve: perturbation result; Dotted curve: exact diagonaliza- 
tion result; Inset: Schematic of the potential along the verti- 
cal (z) direction, (b) Upper panel: Weighted scattering rates 
Ti-,j between different energy levels (from "spin-up" to "spin- 
down") vs. the electric field, r^otal * s ^ e total weighted 
scattering rate from the "spin-up" to the "spin-down" states. 
Lower panel: Energy level ex of the z direction of the double 
QD vs. the electric field. 



one obtains the wave function: 
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in which Ai and Bi are the Airy functions. The coeffi- 
cients together with the eigenenergy e A can be obtained 
from the boundary conditions ip z \(z — ±(a/2 + d)) = 0, 
the continuity conditions at z = ±^a and the condi- 
tion of normalization J ip* x (z)tp z x(z)dz = 1. The elec- 
tron Hamiltonian in the x-y plane is H e = Ho + H so , 
where H = (P 2 + P 2 )/(2ra*) + V c (r) + H B is electron 
Hamiltonian without the spin-orbit interaction, in which 
P = (P x ,P y ) = -ihV + (e/c)A with A = (B/2)(-y,x) 
is the electron momentum operator, to* is the elec- 
tron effective mass. Hb = \g^sBa z is the Zeeman 
energy with a z denoting the Pauli matrix. H so = 
jfe ~}2\{Pz)\{—PxPx+PyOy) is the Dresselhaus spin-orbit 
coupling 22 with (P 2 ) A ee -h 2 Jr zX (z)d 2 /dz 2 ij zX (z)dz 
and 7 = 27.5 A 3 -eV 25 . For the small applied gate 
voltage, the Rashba spin-orbit coupling 23 is unimpor- 
tant in this study 26 . The eigenenergy of H i s E nl(J 
Ml{2n + |/| + 1) — hlu>B + ctEb, in which Q = + lj^ 
lob = eB/(2m*) and E B = ^gusB. The eigenfunc- 

tion (r\nla) = N n ^{ar)\ l \ e -^ 2 ' L l n l {(ar) 2 )e ue X v with 
JV„,I = (a 2 nl/n(n+ \l\)l)^ 2 and a = yJm*Vl/h. is 
the generalized Laguerre polynomial. \<y represents the 
eigenfunction of o z . In these equations n — 0, 1, 2, • • • , 
I = 0, ±1,±2, ••• and u = ±1 are quantum numbers. 
From the eigenfunction of H 0} one can construct the 
wave function |^) of H e by either the perturbation 
calculations 12 ' 14 modified by the right energy corrections 
pointed out by Cheng et al. 15 or the exact diagonalization 
approach. 15 

The SRT r is calculated from t _1 = ^2nf fi^i^f in 
which fi — C exp[— £i/(fcsT)] denotes the Maxwell dis- 
tribution of the i-th level with C standing for the nor- 
malization parameter and 

qAi 

- fiw qAl ) + (n qAl + \)8{E f -E t + fiw qAl )l (3) 



is the transition rate from the i-th level to the /-th one 
due to the electron-phonon scattering due to the defor- 
mation potential with |M qs ;| 2 = KB 2 q/2Dv s i and the 
piezoelectric coupling for the longitudinal phonon mode 
with |M qpi | 2 = {32fnr 2 e 2 e 2 J K 2 Dv sl )[(3q x q y q z ) 2 /q 7 } 
and for the two transverse phonon modes with 



(32h7T 2 e 2 e 2 J K 2 Dv st q 5 )[q 2 x q 2 y + q 2 y q 2 z + 

QzQx ~ (3<fa<Zy9z) 2 /q 2 }- "qAi represents the Bose distribu- 
tion of phonon with mode Ai and momentum q at the 
temperature T. Here S = 7 eV stands for the acous- 
tic deformation potential; D = 5.3 x 10 3 kg/m 3 is the 
GaAs volume density; = 1.41 x 10 9 V/m is the piezo- 
electric constant and k = 12.9 denotes the static dielec- 
tric constant. The acoustic phonon spectra are given by 
w qq; = v siq for the longitudinal mode and ujqpt — v st q 
for the transverse modes with v s i = 5.29 x 10 3 m/s and 
v st = 2.48 x 10 3 m/s being the corresponding sound ve- 
locities. 



The states i and / in Eq. (3) are the eigenstates of the 
Hamiltonian H e . In order to demonstrate the physics 
clearly, we first use the corrected perturbation method 
by Cheng et al. 15 to study the SRT. For the double dot 
system, we need to include the lowest two energy levels of 
z direction which we label as \l z ) and \2 Z ) [Eq. (2)]. In x- 
y plane, the lowest six energy levels of Ho for each QD arc 
considered, i.e., |00+), |00— >, |01+), |01— >, |0-1+), and 
— 1 — }. The wave functions of the lowest four states of 

p2 

H e + ^-7 + V z constructed from these levels are therefore 
given by 
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with the corresponding energies being: 
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In these equations £?„; CTjA = E n ia + e\: B\ = icry A (l - 
eB/(2Ha 2 ))/(E Q -i-.x - £ o+.a) and Ax = ^ 7 *(1 + 
eB/(2Ha 2 ))/(E 01+ ,x - P o-,a) with 7 * = i{P 2 )x/h 2 . 
A(= 1,2) is the quantum number of z-axis. Now we 
calculate the spin-flip rates from the "spin-up" states 
\^2m-i) to the "spin-down" ones \^2m) ( m = 1, 2) due to 
the electron-phonon scattering. There are nine spin-flip 
scattering rates. The scatting rate from the "spin-up" 
state i to the "spin-down" one / reads 

,7T/2 

T^ f = \A f - £,| 2 K + [1 + s S n (* - /)]/2}9 3 / d9 

Jo 

x C LD q 2 sin 3 6 + C LP q 2 sin 7 6 cos 2 6 + C T p sin 5 6 
x (sin 4 9 + 8 cos 4 0)1 e- q2 sin2 e ' 2 1% (q cos 9) \ 2 , (12) 

in which I if (q z ) = {tp z i\e lq!lZ \ip z f) and q = \E t - 
E f \/(hv x a). C LD = E 2 a 3 /(8nhv 2 sl D), C LP = 
9e 2 e 2 4 air/(fiK 2 Dv 2 l ) and Ctp = Tte 2 e\ A a/(hn 2 Dv 2 st ) in 
Eq. (12) are the coefficients from the electron-phonon 
scattering due to the deformation potential and due to 
the piezoelectric coupling for the longitudinal phonon 
mode and two transverse phonon modes respectively. 

In Fig. 1 we plot the SRT of a typical double dot with 
d = 20 nm, a = 10 nm, d = 5 nm, Vq = 0.4 eV and 
B = 0.1 T at T = 4 K as a function of electric field E. 
The solid curve in Fig. 1(a) is the result from the per- 
turbation approach. It is interesting to see that the SRT 
is increased about seven orders of magnitude when the 
electric field is tuned from 0.1 kV/cm to 1.3 kV/cm. The 
physics of such gate- voltage-induced dramatic change can 
be understood as follows: When the gate voltage is small, 
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due to the large well height Vo and/or large inter-dot dis- 
tance a, the electron wavefunction (along the z-axis) of 
the lowest subband of each well is mostly localized in that 
well due to the high barrier between them and hence the 
difference of the lowest two energy levels is very small 
(about eV). When a gate voltage is high enough, 

electron can tunnel through the barrier and the wave- 
functions in the two wells get large overlap. Therefore 
the separation between the lowest two levels e\ and £2 
increases. This can be seen from Fig. 1(b) where the 
energies of the lowest two levels along the z-axis E\ and 
£2 are plotted against electric field E. From Eqs. (8-11) 
one can see that the first two levels (E\ and E2) and 
the next two levels (E3 and E4) are mainly separated 
by the energy along the z-axis, i.e., E\ and £2. Such 
an increase makes the electron-phonon scattering more 
efficient when the energy difference £2 — £1 is not too 
big. Therefore, by applying the gate voltage, one finds 
the SRT first decreases. Nevertheless, with the further 
increase of the gate voltage, half of the lowest four lev- 
els are quickly removed from the spin relaxation chan- 
nel and the SRT is enhanced. As a result, there is a 
minimum of SRT with the gate voltage. This can be 
seen from the same figure where the weighted scattering 
rates = fiT^f) between different levels are plot- 

ted versus the electric field. The leading contribution to 
the total scattering rate comes from f 3 ^ 2 at small field 
regime. When the electric field increases from 0.5 kV/cm 
to 1.3 kV/cm, f 3^2 decreases rapidly due to the separa- 
tion of £\ with the electric field but f i^ 2 keeps almost 
unchanged as both levels E\ and E2 correspond to the 
same lowest level £\ along the z-axis. Finally for large 
field, f i_2 defines the total scattering rate. It is further 
noted that although we performed the average of the ini- 
tial and the sum of the final states in calculating the SRT, 
the leading contribution comes from the scattering from 
E3 to E 2 at low electric field and the scattering from E\ 
to E 2 at large one. 

The large variation of f3_»2 around 1 kV/cm can be 
estimated as following: As the electron-phonon scatter- 
ing due to the piezoelectric coupling of the two trans- 
verse phonon modes is at least one order of magni- 
tude larger than the other modes, we only consider 
the third term in Eq. (12). From our calculation, 
£1 = (3.25 x 10- 4 £/(kV/cm) + 0.15129) eV and £ 2 = 
(1.68 x 10- 3 £/(kV/cm) + 0.1513) eV. The energy split- 
ting between E 2 and E 3 can be approximated by £2 — £i- 
Therefore AE 23 = (1.36 x 10- 3 £/(kV/cm) + 5 x 10" 5 ) 
eV approximately and q = AE 23 / [hv st a) . As the vari- 
ation of |ii2(<;cos6>)| in Eq. (12) is within one order of 
magnitude, we approximately bring it out of the inte- 
gral. Then the remaining integral J^ 2 d9 sin 5 #(sin 4 9 + 

8 cos 4 6)e~ q2 sm2 B l 2 can be carried out analytically: 
\B{\; 5)*(5; -q 2 /2) + 4B(|, 3)$(3; f ; -q 2 /2) with 
B{n;v) and $(a;7;z) being the Beta function and the 
degenerate Hypergeometric function separately. When 
E = 0.1 kV/cm, the value of the integral is 10 _1 and 
when E = 1.3 kV/cm, it becomes 10~ 6 . Meanwhile, 



with the change of the electric field from 0.1 kV/cm to 
1.3 kV/cm, although g 3 |-4/ — £>i| 2 is increased by one 
order of magnitude, I-/23I 2 is decreased by one order of 
magnitude and the distribution function / 3 is decreased 
by another two orders of magnitude. Therefore, f 3 ^ 2 de- 
creases about seven orders of magnitude when E is tuned 
from 0.1 kV/cm to 1.3 kV/cm. 

As pointed out by Cheng et al. 15 and confirmed by 
Destcfani and Ulloa 27 that due to the strong spin-orbit 
coupling, the perturbation approach is inadequate in de- 
scribing the SRT even when the second-order energy cor- 
rections are included. Therefore, in Fig. 1(a) we further 
plot the SRT calculated from the exact diagonalization 
as dotted curve. Similar results are obtained although 
again the SRT from the exact diagonalization approach 
differs from the perturbation one. 
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FIG. 2: SRT calculated from the exact diagonalization ap- 
proach vs. the electric field at (a) different magnetic fields 
with do = 20 nm and (b) QD diameters with B = 0.1 T. In 
the calculation a — 10 nm, d = 5 nm, do — 20 nm, Vo = 0.4 
eV and T = 4 K. 

Now we investigate the magnetic field and dot size de- 
pendence of the SRT in Fig. 2(a) and (b) by exact di- 
agonalization approach. Again one observes a dramatic 
increase of the SRT by tuning the electric field up to 
a certain value and then the SRT is insensitive to the 
electric field. For small dot size (do = 10 nm), one 
even observes a twelve orders of magnitude change of 
the SRT by tuning the gate electric field to 2.6 kV/cm. 
The dramatic variation of the SRT has been explained 
above. Now we discuss why the SRT decreases with 
magnetic field and dot size observed in Fig. 2 in the 
electric-field-insensitive part. From Fig. 1(b) one finds 
Ti^2 is the leading contribution to the total scattering 
rate in this part. The energy splitting between the first 
and the second levels AE12 oc B. As AE12 is about 
10~ 5 eV, n q ~ k B T/AE 12 and n q q 3 oc (AE 12 ) 2 . More- 
over \Ai ~B X \ 2 = {a^AEBUJBf /ih^uil) 2 oc B 4 prox- 
imately. As a result, the coefficient before the integral 
of the electron-transverse phonon scattering due to the 
piezoelectric coupling is proportion to B e . Although the 
integral has a marginal decrease with B, Pi_>2 still in- 
creases with B. Similarly, one can explain the change of 
the SRT with the dot diameter d - 

It is noted that in order to obtain the large variation 
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FIG. 3: SRT calculated from the exact diagonalization ap- 
proach vs. the electric field at (a) different barrier heights Vo 
with the barrier width a = fO nm and (b) different barrier 
widthes a with Vo = 0.4 nm. In the calculation, d — 5 nm, 
do = 20 nm and B = O.f T. T = 4 K. 



of the SRT by a gate voltage, it is important that the 
barrier between the QDs should be large enough so that 
without a gate voltage, the two dots are decoupled (and 
there is no energy splitting along the z-axis) . This can be 
clearly seen from Fig. 3: With the decrease of the barrier 
height Vo or the inter-dot distance a, the tunability of 
the SRT by the gate voltage decreases. 
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FIG. 4: SRT calculated from the exact diagonalization ap- 
proach vs. the electric field at different inter-dot distance a 
with low barrier height Vo = 0.05 eV. In the calculation, d — 5 
nm, d = 20 nm, B = 0.1 T and T = 4 K. 

The double dot system proposed in our scheme can 



be easily realized with the current technology. 28 ' 29 Nev- 
ertheless, it is not essential to use such a high barrier 
height system to obtain the large spin manipulation. For 
ordinary barrier height widely used in the experiment 
(which is about one order of magnitude lower than Vo 
used above), one can still achieve the similar manipula- 
tion by increasing the distance a between the two QDs 
as shown in Fig. 4 where the barrier height Vo = 0.05 
eV. One finds that for small Vo, if the barrier width d is 
large enough, one can still get the large change of SRT. 
Especially, in the case of a — 30 nm, eleven orders of 
magnitude change of SRT is obtained by a small gate 
field. 

In conclusion, we have proposed a feasible scheme to 
manipulate the spin relaxation in GaAs vertical double 
DQs by a small gate voltage. The SRT calculated can 
be tuned up to twelve orders of magnitude by an elec- 
tric field from the gate voltage less than 3 kV/cm. This 
provides a unique way to control the spin relaxation and 
to make spin-based logical gates. The conditions to re- 
alize such a large tunability are addressed. The double 
dot system proposed in our scheme can be easily real- 
ized in the experiment. Finally the proposed large or- 
ders of magnitude change due to the gate voltage will 
not be reduced by the hyperfine interaction with nuclear 
spins 30 ' 31 as the SRT due to this mechanism in our case is 
around 10 3 s at 0.1 T. Finally we point out that differing 
from the earlier reports 32,33 where a strong variation of 
the SRT is obtained from the anticrossing of the energy 
levels induced by the Rashba spin-orbit coupling by in- 
creasing the magnetic field 32 or the inter-dot distance, 33 
there is no anticrossing/crossing of the energy levels in 
our scheme. Moreover, the tunability of the scheme pro- 
posed in the present paper is better as one only need to 
tune a very small gate voltage (to tune the electric field 
from 0.1 to 1.2 kV/cm) to obtain a surge of the SRT up 
to twelve orders of magnitude in contrast to the large 
magnetic field of several tesla to obtain the variation up 
to seven orders of magnitude. 32 
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